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ELECTRO-OPTICAL HSTHOD A?fD SYSTEM FOR DETERMINTVG 
THE DIRSCTIOM OF MOTIOH IN DOU 3 L E--EXPOSURE 
VELOCIMETRY 3Y SHI?TC.VG AM OPTICAL IMAGE FIELD 

The invenrtioa herein resulted from work 
5 which -is supported by TSI Inc. and the National 
Science Foundation under NSF ATM 8^-00509. 

BACKGROUJTD OF THE INyENTIOii • . 

The present invention is an improvemeat on 
that previously described i n • c o- p e.n d ing 
10 application 023,773, of which the present 
application is a can tinua t ion-in-p ar t . The 
benefit of the filing date of that co-pending 
application with respect to coianon su^'ect matter 
is herewith claimed. 

■15 The present invention p e r t a i n s - 1 o. the 

neasureoent .of the velocity of -fluid flows ■a.nd, 
'more particularly, to systens in whi.ch the fluid 
flow rate is inferred- froni neasureaents mads in 
the context of double— pulsed particle . iaage 

. 20 velocioetry , in which saall scattering particles 
are illuatinated by two short pulses of laser 
light or. other light, and their 4jiages are 
record'ed photographlcalLy. th produce a -record 
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from which the particle velocity can be 
. determined by measuring the displacement of the 
particle iaages. 

F.or specific ■ background material,, reference 
5 - say be made to U.S. patent 4,729, 109.. assigned to . 
the assignee- of the present invention; also to 
'~ .co-pending application serial number 023,773, 
..'"also assigned ■ to. the assi'gnee of the present 
invention. 

. 10 ■ As will, be appreciated by the' references 
cited above, particle image velocimetrv fPIV) in 
a well-established technique in experimental 
photo-mechanics .for quantitatively measuring 
velocity data at a given instajit of time over an 
15 extended flow field. A variety of other 
references which wer e • cited*; in co-pending 
•'..application 0^3^, 773- may ' be coft-sulted for the 
" . details of -the 'technique involved. 

Generally speaking, in the PIV technique the 
20 flow is seede'd with smal'l particl'es,' 'typically 
ten micrometers or less in size, and illuminated 
by a thin sheet of pulsed laser light, typically 
coded with a .double pulse. 'Particles moving 
within the light sheet are recorded 
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photographically as pairs of particle iaages. 
The local fluid velocity is found on a grid of 
snail "interrogation spots". on the photograph by 
sequentially measuring the displacement -^^of 
the images vlthin each jsppt. Interrogation is 
accomplished fay a variety of methods including 
the Young's fringe method, direct measurement of 
image-to*- image displacement, and spatial 
correlation analysis. By successively analyzing 
the photograph at many adjacent ^interrogation 
regions, local velocity may be inferred in a 
large nusnber of points lu ^ae flow field. 

Given identical conditions in viewing each 
of the two exposures on a double-pulsed PIV 
photograph, there exist no characteristics on the 
photograph to distinguish first, images from 
-second isages. As a result, measurement .of the 
particle image separation cannot detera'ine the 
polarity of the fluid velocity and the velocity 
v^ldtor is ambiguous in sign: ±\i^ Clearly, it is 
necessary to determine the order of the tvo 
exposures on .the photograph to e.l lain ate 
directional ambiguity. 
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- In "accordance with the invention described 
^ in co-pending application 023,773, a techniqua . 
known as "iaase shiftiag" has been shown to-be 
■ effective' in resolving 'directional ambiguity. 
5 This Qethod displaces the photographic iniase 
. field by an appropriate uniform, kno'-vn distance 
between the '"first and second 'pulses .' As a 
result, the second image of each .particle laage 
pair is shifted by a difiplace.oent such that . 
10 the most negative fluid velocity still produces a 
positive displacement o.f the second particle 
isifige with respect to ths first.. Alfter'the 
in-terrogation , the artificial shift 
subtracted laathesiatlcally to -obtain the actirsl 
■15 fluid velocity. 

Several.- methods have bedn ' irapleaented- to - 
produce the -laage shifting during the recording' 
procedure, in particle laage ' velociaetry . laasa 
shifting has been accomplished by placing a 
20" rotating . airro.r- .i.n front of the photographic, 
lens, and by translating the camera apparatus 
between exposures. . Other aethods have- also been- 
proposed, including noviag the-fila within the 
camera -and translating the 'caraera lens with a 
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piezoelectric device. The drawback with each of 
the above-noted methods is_ the mechanical motion 
involved in part of the recording apparatus 
between the first and second lighf pu.lses . 
5 Consequently, they are subject to limitations on 
speed due to the dynamics of moving components, 
and limitations on accuracy and consistency. 

Accordingly, it is a' primary object' of the 
present invention to overcome the limitations on 
10 operations at high and low speeds in connection 
with image shifting in pulsed velocimetry 
systems, and to carry oiit the operations with a 
high degree of accuracy. 

Ano'ther obje'ct is to' permit image ' shifting 
15 in periods .of tine less than one microsecond. 
- - -• A-fMrther • object is to provide precise 
•reproducible shifts without the nee-d for moving 
parts to accomplish them. 

- SUMMARY OF THE INVESTIOM 
20 The above and other objects are fulfilled in 

accordance; with the present invention by a 
primary feature involving an electro-optic iaags 
shifting arrangement - including a biref ringent , 
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uniaxial crystal plat.e; and a camera lens 
interposed be.tveen said crystal plate '■ and 'the 
image plane, the crystal plate being cut.- with its 
. "surface nornalV aligned parallel, to t.he 

'5 crystal's optic axis, s'uch that a given light ray 
from any source* in the object plane' is refracted 
into tvo parallel light rays having mutaally 
orthogonal states of linear polarization. 
Accordingly, for a. given , scattering particle in 

id .the object plane, tvo iaages of the particle are 
formed>at the film plane of the camera separated 
by k'g, which is the anount of iaage shifting 
introduced by the present invention. Howeve.r, it 
will be understood that when the moving particle 

15 or the iaage thereof is recorded by. means of two 
illuminating pulses, four iaages are thereby 
formed'. ' ;*' 

In accordance with a more specific feature 
pf the present invention,, means are provided for 

20 controlling the polarization of the' light" 
scattered from the seeding particles so that the 
first illuminating pulse from the source scatters 
light which is linearly polarized and parallel to 
the principal plane of the birefringent crystal, 
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and the second pulse scatters light which is 
linearly polarized .at 90 degrees, with respect to 
■ the light from the first pulse, ^whereby only two 
particle iaages will be -'forned for_ a single 
.5 moving particle. 

Other and further obj-ects, advantages and 
features of.. tHe^ present invention will be 
understood by reference to • the -following 
specification in conjunction with the annexed 
^0 drawing,' wherein like parts have been given like. ^ 
nuabers. 
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" BRIEF DESCRIPTION OF DRAWIXG " . • - 

Fig. 1 -is a schematic diagraa of the 
. electro-optic image shifting system of the 
- .present invention. . 

5 F.ig. 2 ilriustrates "the polarization of 

linearly-polarized light at 90 degrees by a 
spherical particle. ' - _ • 

-Figs. 3A and 3B ilLustrate the light 
scattering response of plastic particles inwatpr 
10 (n/nosl.20) for linearly-polarized illuoination 
694.3nm),'at 90 degrees si.de-scatter , using 
a finite light gatherino mne- Kuaeri-cal 
calculations incorporate Mie scattering theory 
and assume a lens aperture of fl6. . The .two 
15 horizontal bars' -in each graph bound the 
■approximate range of scattering response suitable 
for Kodak Technicar Pan 2415 film. . . 

. F.igs.'4.4, 4B, and 4'C.. ilia's trate image shift 
response using a calcite crystal (hq » 1.65259, 
20 ne = 1..4S382) , having an optic" axis angle ^3 = 
53.53 . Open circles indicate observed - 
measurements. Fig. 4A shovs principal iaage shift 
vs. the tangent of the pr in c i p a 1 an g u 1 ar 
derivation from the photographic axis; FigTAS is 
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the sane as Fig. 4A tut reduced . scale; Fig. 4C 
shows transverse image shift vs. the tangent of 
the transverse angular, deviation from the 
photographic axis. -Xgo is the on-a:tis image- 
shift, and Si is the iaage distance. . 

Figs. 3A and 5B illustrate PIV velocity 
vector' maps of a. unif oraily-displaced water- flow 
field "Slaving a s'nall transient flow velocity," 
using eleictr.o-optici isiage shifting. Fig,"5.\ 
shows stationary reference frame. _ Fig. 5B shows 
reference frane noving at the uniform translation 
velocity, 

DESCRIPTION OF PREFERRED EMBODI.MENT 

Referring now to the figures of the drawing 
15 and, in particular. Fig. 1, there is illustrated 
in .this figure an espbodiment 'of the electro-optic 
image shifting system of the present invention. 
As will be seen therein, a birefringent , uniaxial 
crystal plate 10 ( f or -exam.p le , calcite) is 
■ 20 - positioned, in .front o-f a camera lens i2. , The 
. plate is cut with its "surface normal" aligned 
-parallel to the crystal's principal plane and 'at 
an obliqiie angle to the crystal's optic axis, 
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such that- a given light. ray from any source in 
the object plane 14 is refracted into two 
parallel light rays 16a and 16B, coaaonly 
referred to as an b-ray and an -e-ray* having 
S.Butjially orthogonal states of linear' 
polajlzation. Consequently, the separation of 
the" tvo rays is nearly uniform for a wide range 
of incidenX angles and scaled to the desired 
image shift Xg seen on the image plane 13 for the 
10 chosen camera magnif icatj.on. ' 

For-' a given scattering particle 20 in the 
'. . ooject. plane 14 within the flow fxelfl, tvD i:ra3es 

of that particle are foraed at the filn or iaage 
plane- IS Qf the cajnera, separated by Xg. The 
15 Intensities of the two images are equal to the 
* _ ■• intensities of the two orthogonal components of 

■ -• polarization, with respect to the crystal's 

principal, plane, of the light incident' on the 
plate 10. Accordingly, when a moving particle is 
■- - ■•-20 recorded with two illuniaating pulses, f our _ 
iaages ere formed, arranged in 'a parallelogram 
ha.ving two sides equal in length to the inage 
displacement [AjC \ and two sides equal to\X^ . 
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If the • polarization of ' the light scattered 
from the seeding particles is controlled ' so that 
the first illuminating pulse scatters light which 
is linearlj " polarized and parallel to the 
5 principal plane of the birefringeat crystal plate 
10, and the second pulse gathers light which is 
linearly polarized- at -90 degrees with respect to 
the light from the first pulse,' then only two 
particle images will be formed for a single 
10 particle, displaced by the vector sun of Ax and 
JCs_ . This is the desired result for any image 
shifting technique. 

In the system shown in Fig. 1, the 
polarization of fhe scattered light is controlled 
15 by switching the polarization •• of - the laser beaa • 
30 at tTie' exit of the' cavity' of double-pulsed* . 
laser 32. A Pockels cell 34, well-khown per se, • 
is used with the double-pulsed laser to switch 
the .linear polari.zatio'n, between pulses, from a 
'20 vertical state to a horizontal state (or vice - 
versa). The output of Pockels cell 34 is 
transaitCed to sheet forming optics means 36, 
also well-known, thus producing the "light sheet 
38. An alternate 'systea"to that of Fig. 1 
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. incor porafes the output from .two orthogonally- 
polarized single-pulsed lasers, fired in 
sequence, which are coralined using a polarizing • 
beaa-splitti'ng cube. 
5' .In the system of Fig. l; the extinction by;' 

the uniaxial pXate 10 of the shifted image of the 
first exposure and the- unshifted iaage of the 
second exposure depends upon creating linear, 
orthogonal" states of polarization over the solid 
10 angle Jldefined by the photographic aperture. In 
general, for spherical particles, the 
polarization &£ the sc^tc^red wnve >. t . rs.-. . by 
the Mie scattering coefficient 

where A and 3 are cooplex functions of ^and all 
15. other quantiti.es .are de-fined 'in Fig. 2. When ' 
side-scattering pccurs over 'an inf initesiaal 
solid angle, perpendicular to the direction of 
illumination^^ , e-S- centered on the 
photographic axis. z,.the po-lariz-ation of the 
20 scattered light. wave is linear. If the 
illusinating beam is polarized vertically. (Eqi), 
then the side-scattered light is also vertically - 
: polarized: (J^J s^^^c^q')^^ If -the .illuainatiag 
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wave is polarized horizontally, then the side- 
scattered wave is also polarized horizontally: 
s-^WjE.' ^° ideal limit of 

.". " 'infinitesimal solid angle, the first and second 
5 exposures of a particle centered in the camera's 
"field of view would ba polarized, respectively, 
. in the X and Y d i r e c t i o n s . of the camera 
coordinates .(Fig. • 1), as desired.' 

• Practical applications of this approach 
10 require . orthogonal polarization of the two 
scattered light waves when the photographic field 
of' view is finite and the solid angle of the lens 
is finite. The latter.'ef f ect has been evaluated 
by calculating the -scattering" characteristics of 
15 typical PIV seeding particles, using Mie 
scattering theory. - 

■ Figure 3 .shows the results of calculati-ons 
of the scattering efficiency (defined as the 
. ratio of the- particle iaage intensity to the 
.20- illumination -i-n tensity) p_erforaed for a - 
polystyrene particle In water, in which the 
particle diameter is varied u?' to 10 n and a' 
lens aperture of f 16" is used. In' Fig , 3A , • two 
.curves are shown corresponding to the X and Y " ' 
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•coaponents of the scattered light, for a "y- 
polariz'ed light source. The Jf, polarization 
coaponent is over four orders of oagnltude larger ' 
than the X conponent for the entire range of 
5 particle diameters. In Fig. 3B, the source is 
switched '-to z^polar ization , and the ^ and Jf^ 
/'respo-nses likewise switch. The X-coraponent 
' doainate.s .by over four .orders of magnitude. 
-. Calculations using a variety of other particle 
10 sizes and coapositions (including hollow spheres 
and metallic-coated spheres) show similar 
tendencies for the scattere'd light waves to 
retain linear, orthrogonal polarization when the 
scattering occurs over finite solid aagles normal 
15 to the illusinating light. 

The effect of finite fitsld- 0;"f .view is ■ to 
include- scattering from off-axis particles whose 
mean scattering angles are not perpendicular to 
the light sheet (e.g.-fl h ?Cf ) . In general, the 
- 20. polapizations at. th^se angles need not be linear 
or orthogonal. However, experimental -exaainatioa 
of the ratio of the horizontally "and vertically 
polarized 'waves indicates that the extinction 
. ratio is large over.'a i20* f ield-of-view for 13 
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aicron polystyrene particles in -water and 4 
aicron silver-coated glass spheres in water. The 
use of a finite ' field of view is thus not a 
limiting factor in these cases. 
5 The design of the itniaxial crystal plate 10 

(Fig. 1') involves appropriate selection of the 
angle between the- optic axis -0^ .'and the noraal to 
the transaittlitg surfaces ,-- and' the. 'plate 
thickness T. is chosen so that the image 

10 shift (X,^ Y) is nearly uniform across the 
photograph, and T is, selected so that the shift 
is scaled -to the desired- constant value X3 for 
the user-chosen camera magnification. The 
optimal value of. has -been determined 

15 computationally by modeling the propagation of 
light through -uniaxial calcite material using 
Huygetis's vavefront" construction in .three 
dimensions.' A" value » 53.58 degrees results 
in a shi£t magnitude Xg along the principal plane 

20. of' the- cry stal ■ which is. very uniform. The ' 
calculations in Fig. 4A indicate -a maximum 
variation of the shift in the X-direction that is 
less than 0.15 per cent. -The computations 
indicate the presence of a non-negligible 
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• • ■ ■ * 

transverse component Ys(X,Y), perpendicular to 
the .principal plane, which is shown in Fig. 4C. 
The image shift vector ^1^^ is thus " sonewhat non- 
uniform, but. the non-uniformitj is known a priori 
5 and thus may be corrected computationally during 
the interrogation procedure* 

• EXPERIMENTAL RESULTS . . . • - " 

An electro-optic iaagf shifting system using 
• the- Pockels cell system of Fig. 1 has-been 
10 constructed. The response "time of the Poc-kels 
. cell is approximately 10 ns, which is coaparable 
to the pulse width of the double pulsed lasers 
used forPIV, and therefor^ capable of iaage 
shifting the fastest flow fields that are 
-15 currently, contemplated. Currently - there" fire 
three calcite plates- in use having thicknesses of 
2 mm, 6 -aa, and 10 mn, 7 which oay be mouated 
separately -or combined- in front of the.caaera 
-.lens- to- provide-image shifts . magnitudes oi 
20 0.2m3<Xs<1.8 mm. . 

Two- experiments were performed to evaluate' 
•.. the imase-'shif ting technique. "In the first, the 
accuracy of "image-shifted aeasurements vas 
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evaluated across a PIV photograph. For this- 
esperiment, single-palse illuaination was used to 
record a stationary flow field' illuninated by a 
circularly polarized beam. The Pockels cell Has 
5 not activated in this -experiment. As a result, 
the PIV photographic field contained two images 
of each particle, separated- by the Local value of 
J^.^°^- i"*S.e pair. Over three hundred 
particle iaage displaceaents. were neasured across 
10 the . pho t.ograph . -A selection of. these 
laeasureoents ,^ collapsed along the X- and 
directions, are plotted' in Figs. 4B and 4C. The 
aeasureaents agree veil, with the nun.erical 
prejiictions-. The RMS error of the neasureaeats 
15 is approxiaately 0.8% full-scale, most" of vhich 
may be attributed to the -accuracy limitations of • 
.the PIV interrogation syscen ' used . to. -extract 
inf ofaa'tion from the photograph. - 

In the second experiaent, the entire i3age 
.. .2D -shifting .systen was evaluated using PIV photos of- 
nearly-quiesce'nt waxer contained in a test 
■ section- which was horizontally displaced wi-th a 
uniform ..velocity U. A small transient -velocity 
' u^Cx.y.z) was present in the' flow .during- the 
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uniform translation, where \ut\a»2e < 0.15 U» The 
calcite plate was oriented to provide a vertical 
image shift. The photographs were taken over an 
extended field of view (±15°) using, a finite 
5 camera aperture (f8). Each interrogation of a 
PIV photograph of the flow generally produced a 
highly detachabl-e neasurement of the particle 
- image displacement, consisting of the vector sua 
.of the iaage shift JCs(X,Y) , the displacement 
10 littV, and a .small component which is attributed 
to the .transient displacement M^tu is the 
masaification 'of the camera). The vector liap of 
the interrogated PIV photograph (Fig. 5A) shows a 
-'highly resolved- velocity f ield ]u + jit(x , y) . When 
15 the translation velocity ^ is subtracted from the 
vector field, Ithe"' resulting map (F.ig. SB) clearly 
. _ reveals the tr3nsi"ent structure! 

. SUMMARY AND CONCLUSIONS 

In suanary^ "image shifting based on electro- 
20 optic" hardware is an effective approach in PIV 
that offers several advantages over nechanical 
shifting techniques.- Two systeras" for 
' implementing this approach, one iocor porating a 
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. Pockels cell to modulate the output froo a single 
. laser and the other incorporating the output f roa 
two orthogonally polarized lasers have been 
described". - -The use of a properly designed 
5 calcite eleoent results in an iaage shift which 
is nearly uniforra over a vide light-gathering 
.' cone. The aethod is accurate, reliable and fast 
■ enough to a.c c om nio d a t e the h i g h e s t • f 1 u i d' 
- velocities likely to be studied with particle 
10 image velocinetry , and it is equally applicable 
to laser speckle velociaetry. 

It is understood that the effectiveness of 
either of the two systeias described above depends 
upon retaining- light polarrization af ter 
13 sidescatter. Consequently they-may be limited to 
-certain ranges of particle types,- particle .sizes, 
■angular fields of view and caaera apertures-. 
Further investigation is required to determine 
these limitations. An alternative electro- 
20 .optical system, vh^ch incorporates a wide-angle 
Pockels cell in front' of the uniaxial plate for 
directly acdulating the particle-scattered light, 
is anticipated . to" expand the effectiveness of ■ 



SUBSTITUTE SHEET 



wo 90/10876 



PCr/US90/0I322 



20 

electro-optic- iraage shifting, and is curr.eatly 
being' investigated. 

While there has been shovn and described 
vhat is considered at present to be the preferred. 
5 esibodiment of the present invention, it -will be 
appreciated by those skilled in the art that 
■modifications of such- eabodiaent nay be made. It 
is therefore desired, that the invention not be 
limited to this embodiment, and it is intended to 
10 cover in the appended claims all such 
modifications as fall within the true spirit and 
scope of the invention. 
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CLAIMS 

.1. In a particle -image velocimetry system, 
in which images are created by accurately tioed 
puls-es of .light f r on - a li gh t source', and 
successive particle* images, resulting froa' 
seeding .particles being displaced by the motion' 
. o£ a fluid, are . recorded ,. such that the fluid- 
velocity is inferred from. the. jdisplacement of the 
image field between exposures ,• the iaproveaenc 
which comprises: 

means' for determining unambiguously the • 
direction of displacement and hence the sign of 
the velocity vector, said means including means 
for shifting the successive particle images vith',', 
r,es>&ct"^o the..f irist image so that all image, 
displacements are positive ' even though physl-cal 
displacements 'o.f the seeding particles may be • 
negative,- said means for shifting being an 
electro-optic image shifting arrangement, 
including a biref ri-ngent i uniaxial .crystal plate;- 
and a camera lens . interposed between. said plate 
and the image planes 
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2. The systea as defined in claim 1, in 
which said crystal plate is cut with- its surface . 
normal aligned parallel to the crystal's 
principal plane' and at an oblique angle "0^ to 
the crystal's optic axis, such that a given light 
ray' from, any source in the optic plane is 
refracted, into two par.allel light rays haying 
a u t li a 1 i y "orthogonal, s t a t e'.s o f . . "1 i n e a r 
polarization. - . 

* 3. In a system as defined in claim 1 in 
which the light source is a double-pulsed -laser^ ^^■''■"•j 
and further comprising a Pockels cell7 sheet- ■ /. 
forming optics^ and a light sheet produced by . 
said sheet-forming optics. / - 

4. In a system as defined .in claim 3t the 
improvement further comprising means for 
controlling the polarization -of _the light 
scattered from the seeding -particles so that .the 
first illuminating pulse scatters light which. is 
linearly polarized and parallel to the principal 
plane "of the birefringent crystal p.late, the 
second pulse ' scatters light which is linearly 
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polarized at 90 degrees with respect to light 

from the first pulse, whereby only • two particle 
-images are formed" for a single particle, 

displaced by the vector sum ofA.x and Xg, where 
Ax is the displacement between successive .images 

and A and x"a is the image shift for each of "the 

images. " 

"5- In a system as defined in claim 4, in 
which said means for controlling the polarization 
of the scattered light includes means for ■ 
switching the polarization of the laser, beam at 
the exit of the laser cavity. 

• 6. . In a system as defined in claim 5, in 
"which said means for switching 'comprises sz 

Pockels" cell "for switching the'l'inear • 
.polarization, between pulses, from a vertical 

state to a. horizontal state. 



SUBSTITUTE SHEET 



wo 90/10876 



1/4 



PCr/US90/01322 




SUBSTITUTE SHEET 



wo 90/10876 



2/4 



PCr/US90/01322 




SUBSTITUTE SHEET 



wo 90/10876 



3/4 



PCr/US90/01322 




SUBSTITUTE SHEET 



wo 90/10876 



4/4 



PCr/US90/01322 




SUBSTITUTE SHEET 



INTERNATIONAL SEARCH REPORT 

CLASSinCATlOW OF SUBJgCT WWTTgW lit lewral gi«...<i>...... . .. 



Aceofomg to im.rnai.oral Pa(ent ClM«fic*t.on (IPC) or to Boin NaKonal Cla»|.fic»t.on ar 

IPC(5) GOIP 3/36 
U.S. CL. 356/28; 73/861.05 



Minimum Docjmeniiiiion Searcftea ' 



Classiiication Syrioois 



356/28, 28.5 
73/861.05, 861.06 



III. DOCUMBMTS COIIHOmiD TO M miJVAWT « 





C tatfon ol Doeumonl. <* wiili indi 


^lion, whore aeoraarialo, of tn* rtlt«tnt patiagts 


RalovanttoCtiimNo 


A,P 


US, A, 4,851,697 SCHO 


DL published 25 July 1989 


1-6 


A 


US, A, A, 733,962 BBEN 


DMDEHL piibHshed 29 March 1988 


1-6 




IV. ciRTincATieM 


Data e( tho Actual Complallan el Ui* imamaMnal S 

21 m 1990 


■weh Data e< MaiHng of tliit InlafliatioiMl Soareh Raoert 


InMfnalional Saarcnmg Authenly 

ISA /as 


C. BUCZINSKI 



